Exclusive Production of the X(3872) in B Meson Decay 



Eric Braaten and Masaoki Kusunoki 

Physics Department, Ohio State University, Columbus, Ohio 43210, USA 

(Dated: February 2, 2008) 

Abstract 

If the recently-discovered ciiarmonium-like state X(3872) is a loosely-bound S-wave molecule of 
the charm mesons D^D*^ or D*^D^, it can be produced through the weak decay of the B meson 
into D^D*^K or D*^D^K followed by the coalescence of the charm mesons at a long-distance scale 
set by the scattering length of the charm mesons. The long-distance factors in the amplitude for 
the decay B — > XK are determined by the binding energy of X, while the short-distance factors are 
essentially determined by the amplitudes for B —>■ D^D*^K and B — > D*^D^K near the thresholds 
for the charm mesons. We obtain a crude determination of the short-distance amplitudes by 
analyzing data from the Babar collaboration on the branching fractions for B D^*^D^*^K using 
a factorization assumption, heavy quark symmetry, and isospin symmetry. The resulting order- 
of-magnitude estimate of the branching fraction for B~^ XK^ is compatible with observations 
provided that J/V' vr+vr" is a major decay mode of the X. The branching fraction for B^ XK^ is 
predicted to be suppressed by more than an order of magnitude compared to that for B^ XK^ . 

PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a narrow charmonium-like resonance near 3872 MeV discovered by the 
Belle collaboration in electron-positron coUisions through the S-meson decay B"^ — > XK"^ 
followed by the decay X — > J/ijj7i~^7r~ [1]. This state has been confirmed by the CDF [2] 
and DO [3] collaborations through its inclusive production in proton-antiproton collisions. 
The discovery mode XK^ has also been confirmed by the Babar collaboration [4]. 

The combined measurement of the mass of the X is 3871.9 ± 0.5 MeV [5]. The width of the 
X is less than 2.3 MeV at 90% C.L. [1], which is much narrower than other charmonium 
states above the DD threshold. The product of the branching fractions associated with the 
discovery channel has been measured by the Belle and Babar collaborations: [1, 4, 5] 

Br[B+ XK+] Br[X J/i/j n+n-] = (1.3 ± 0.3) x 10"^ (1) 

The Belle collaboration recently observed the X(3872) in a second decay mode: X —>■ 
J/'07r+7r~7r° [6]. The invariant mass distribution of the three pions is dominated by a 
virtual cu resonance. The branching ratio relative to the discovery decay channel is [6] 



Upper limits have been placed on the branching fractions for other decay modes of the X, 
including D°D°, D+D', D^D^t:^ [7], Xci7, Xc27, ^/^7, J/^7r°7r° [6], and J/^pr] [8]. Upper 
hmits have also been placed on the partial widths for the decay of X(3872) into e+e~ [9, 10] 
and into 77 [10]. 

The most plausible interpretations of the X(3872) are a charmonium state with con- 
stituents cc [11-13] or a hadronic molecule with constituents DD* [14-19]. Other proposed 
interpretations include an S-wave threshold enhancement in D^D*^ scattering [20], a "hybrid 
charmonium" state with constituents ccg [21], a vector glueball with a small admixture of 
charmonium states [22], and a diquark-antidiquark bound state with constituents cucu [23]. 
Measurements of the decays of the X can be used to determine its quantum numbers and 
narrow down these options [24-27]. The charmonium options include members of the mul- 
tiplet of the first radial excitation of P-wave charmonium and the multiplet of ground-state 
D-wave charmonium. The hc{2P), Xco(2P), and ipillD), whose respective J^*" quantum 
numbers are 1"^", O^"*", and 1 , have already been ruled out as candidates for X. Evidence 
disfavoring each of the remaining states has been steadily accumulating [5, 13]. 

The option of a DD* molecule is motivated by the proximity of the X to the threshold 
for D^D*^: moo+mD^o = 3871.3 ± 1.0 MeV.^ If X is a D^D*^/D*^D^ molecule, its binding 
energy Ei, = rrij^o + m£)*o — mx is = —0.6 ±1.1 MeV. The central value corresponds to a 
resonance, but the error bar allows a bound state with Eh > 0. The possibility that charm 
mesons might form molecular states was first considered some time ago [28-31]. Tornqvist 
pointed out that DD* molecules were likely to be very loosely bound [32]. If the binding 
of D^D*^ or D*'^D^ is due to pion exchange, the most favorable channels arc S-wavc with 
quantum numbers J^'~^ = and P-wave with quantum numbers O^"*" [14, 15, 32]. In 
a model by Swanson that includes pion exchange at long distances and quark exchange 



^ The uncertainty in m^j.o — m^jo is only 0.07 MeV, so the uncertainty in the threshold energy comes almost 
entirely from the 0.5 MeV uncertainty in m/jo. 
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at short distances, the only bound state is in the l"'"^ channel [19]. Another mechanism 
for generating a DD* molecule is the accidental fine-tuning of the mass of the hc{2P) or 
Xci(2P) to the DD* threshold which creates a molecule with quantum numbers 1^ or l^^, 
respectively [18]. 

In the decay X — > J/'07r''"7r~7r°, the invariant mass distribution of the three pions is 
dominated by a virtual uj resonance. In the decay X — * J/ip 7i~^7r~ , the invariant mass dis- 
tribution of the two pions seems to peak near the upper endpoint [6], which suggests that 
they come from the decay of a virtual p° resonance. If the observed decay modes have been 
correctly interpreted as X ^ J/t/jp* and X — > J/'^uj*, the roughly equal branching frac- 
tions in Eq. (2) implies large isospin violations. This immediately rules out all charmonium 
options with the exceptions of hc{2P) and Xci(2P). These charmonium states are excep- 
tional, because if either of their masses was somehow tuned sufficiently close to the D^D*^ 
threshold, the resonant interactions with S-wave DD* scattering states would transform the 
charmonium state into a state that is predominantly a DD* molecule. If the X(3872) is a 
loosely-bound DD* molecule, large isospin violations are expected. They arise simply from 
the fact that the mass of the X is much closer to the D^D*^ threshold than to the D^D*~ 
threshold at 3879.4 ± 1.0 MeV. 

The interpretation of X(3872) as a loosely-bound DD* molecule has important imph- 
cations for the production of the X(3872) [33-36]. If the X is a loosely-bound S-wave 
D^D*^/D*^D^ molecule, it can be produced in any high-energy collision that can create D^ 
and or D*^ and D^. If these charm mesons are produced with sufficiently small rela- 
tive momentum, they can subsequently coalesce into the X. In Ref. [33], this coalescence 
mechanism was applied to the exclusive decay process T(4S') — > Xh^hr . where and h~ 
are light hadrons. This decay can proceed through the decay of T(4S') into a virtual 
and a virtual B~ , followed by the decays — > D^h^ and B~ —>■ D*^h~ or by the decays 
S+ D*^h'^ and B~ — > D^h~, and then finally by the coalescence of the charm mesons into 
the X. Remarkably, the rate for this process can be calculated in terms of hadron masses 
and the width of the B meson only. Unfortunately, it is many orders of magnitude too 
small to ever be observed in an experiment. The coalcscncc mechanism was applied to the 
discovery mode B^ — >■ XK^ in Ref. [34]. An order-of-niagnitude estimate of the branching 
fraction was found to be consistent with current experimental observations if the charge 
conjugation quantum number of the X is C = + and if J/i/j7r~^7r~ is one of its major decay 
modes. It was pointed out that the molecular interpretation of the X can be confirmed by 
the observation of a peak in the D^D*^ invariant mass distribution just above the D^D*^ 
threshold in the decay S+ D^D*^K+. 

In this paper, we present a more thorough treatment of the coalesence mechanism for the 
production of X in the exclusive decays B XK. We extend our previous work to the 
decay mode B^ — > XK^, which has not yet been observed. Our analysis indicates that the 
branching fraction for B^ — > XK^ should be suppressed relative to that for B^ — > XK^ 
by more than an order of magnitude. A measurement or upper bound on the branching 
fraction for B^ XK^ that is much smaller than that for B^ — > XK^ would support the 
identification of the X(3872) as a DD* molecule. 

II. UNIVERSALITY AND THE X(3872) 

The measurement of the mass of the X indicates that its binding energy is likely to 
be within 1 MeV of the D^D*^ threshold. This binding energy is small compared to the 
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natural energy scale associated with pion exchange: m^/(2/i£)D*) ~ 10 MeV, where udd* is 
the reduced mass of the and D*^. The unnaturally small binding energy Ejy implies that 
the S-wave D^D*^ scattering length is large compared to the natural scale l/mT^. We assume 
that the D^D*^ /D*^D^ system has a large scattering length a in the C — -\- channel and 
that the scattering length in the C = — channel is negligible in comparison. The scattering 
lengths for elastic D^D*^ scattering and for elastic D*^D^ scattering are then both 
The large scattering length a could be tuned to +00 by adjusting a short-distance parameter 
in QCD. One possible choice for this parameter is the mass of the up quark, since the precise 
value of the D^D*^ threshold energy depends on this mass. Nonrelativistic few-body systems 
with short-range interactions and a large scattering length have universal properties that 
depend on the scattering length but are otherwise insensitive to details at distances small 
compared to a [37]. Thus, the D^D*^/D*^D^ system will have universal properties that are 
insensitive to any of the shorter distance scales of QCD. 

If the scattering length a is large and positive, the simplest prediction of universality is 
that there is a loosely-bound 2-body bound state. In the case of the D^D*^/D*^D^ system, 
this bound state can be identified with the X(3872). The universal formula for its binding 
energy is 

Eb = 7r^2- (3) 

The reduced mass iidd* is well- approximated by 

mDorriD'O , , 

HDD* ^ ■ (4) 

mx 

The wavefunction of the X for DD* separations r ^ ^/frin is universal. The normalized 
wavefunction is 

^(r) = {2nar/^'-^^tlM. (5) 

r 

Its Fourier transform is 

~, , (87r/a)V2 

The momentum-space wavefunction has this universal form for q <^ m^. 

In the measurement of the binding energy, Ei, = —0.6 ± 1.1 MeV, the central value is 

negative, corresponding to a resonance in D^D*^ or D*^D^ scattering rather than a bound 
state. The small negative binding energy requires a large negative scattering length. If this 
resonance was indeed the X(3872), its dominant decay modes would be Z)*^'/)*" and D*^D^. 
After the subsequent decay of the D*^ or the ultimate final state is D°D°7r° or D°D°^. 
The Belle collaboration has set an upper limit on the product of the branching fractions 
for XK^ and X — * D^D^n^ [7]. Dividing by the central value of the corresponding 

product of branching fractions for the discovery mode, we obtain the limit 

^ ^ <5. (7) 



Bt[X J/i/jTr+TT- 



^ In previous papers [18, 33-35], we also denoted the scattering length in the C — + channel by a, but we 
referred to it incorrectly as the D^D*^ scattering length. 
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This loose upper bound seems to be the only quantitative evidence against the identification 
of X(3872) as a D^D*^/D*^D^ resonance. 

Throughout most of this paper, we will assume that X is a D^D*^ / D*^D^ bound state, 
which requires that a be large and positive. The state of the X can be written schematically 
as 

\X) = {\D'D*') + \D*'D')) + ZT\H). (8) 

where Z is the probability for the X to be in the D^D*^ / D*^D^ state, Zh is the probability 
for the X to be in another hadronic state H, and Z + Y1ih ~ 1- The other hadronic states 
H could include charmonium states such as Xci(2-P), scattering states of charm mesons such 
as D~^D*~ , and scattering states of a charmonium and a light hadron such as J/ip p or J/ipuj. 
Universality implies that Zh scales like 1/a and Z approaches 1 as a increases [18]. In the 
limit a — > oo, the state is a pure D^D*^ / D*^D^ molecule. 

The amplitudes for the scattering of D^D*'^ (or D'^D*'^) with sufficiently small relative 
momentum are universal [18]. If the and D*^ have momenta ±q in the DD* rest frame 
with \q\ <^ ttItt, the universal expressions for their amplitudes to scatter into D^D*^ and 
are 

—1/a — iq 
—1/a — iq 

We have implicitly assumed that the D* has the same polarization vector in the initial 
and final states. The numerator in Eqs. (9) differs from the product of the standard non- 

relativistic factor An/fioD* and the relativistic normalization factor ArriDfnD* by a factor of 
(l/\/2)^ from the projection of the initial and final DD* states onto the C = + channel. 
The amplitudes for D^D*^ or D*^D^ with |q| <^ ttIt^ to coalesce into X are also universal 
[18]: 

A[D'^D*'^^X] = {IGTrmj^/pDD^ay/^, (10a) 
A[D*'^D^^X] = (167rm|//XDD-a)i/2_ ^^qI^) 

We have implicitly assumed that the X has the same polarization vector as the D*. Other- 
wise these amplitudes have an additional factor of e^^ -e^j. . We have simplified the expressions 
for the amplitudes given in Ref. [33] by setting Z — 1 and using the expression for the re- 
duced mass in Eq. (4). 

Any production process for X necessarily involves a wide range of important momentum 
scales. In the decay B — > XK, the largest momentum scale is the mass Mw ~ 80 GeV of 
the W boson that mediates the quark decay process b — > ccs. The next largest momentum 
scale is — 2mc 1.5 GeV, which is the scale of the energy of the s that recoils against 
the cc system. Then there is the scale Aqcd ~ 300 MeV associated with the wavefunctions 
of light quarks in the hadrons and the scale 171,^ ^ 140 MeV associated with the pion- 
exchange interaction between the charm mesons. The smallest momentum scale is set by 
the wavefunction of the DD* molecule: 1/a < 30 MeV if Eb < 0.5 MeV. This hierarchy of 
momentum scales can be summarized by the inequalities 

1/a < < Aqcd -^rrib- 2mc < Mw (11) 
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The two largest momentum scales, whose ratio is about 50, provides a very large hierarchy. 
It can be exploited by replacing the effects of the virtual W by an effective weak hamiltonian 
that involves local interactions between quark fields. The two smallest momentum scales 
in Eq. (11) may also provide a very large hierarchy. The ratio rriT^a is greater than 4.3 if 
El, < 0.5 MeV and it could be much greater. 

We can exploit the universality of nonrelativistic particles with large scattering lengths 
by introducing an ultraviolet cutoff A on the relative momentum q of the or in the 
DD* rest frame. If we choose this cutoff in the range 

l/a<A<m^, (12) 

the behavior of D^D*^ or D*^D'^ with |qf| < A is governed by universality. We will refer to 
processes involving such small momenta as long-distance, while processes involving momenta 
satisfying \q\ > A will be referred to as short- distance. In the next two sections, we apply 
this separation of scales to the decays B — > D'^D*'^K and B — > XK. 

The ultraviolet momentum cutoff A imphes an ultraviolet energy cutoff /{2hdd*)- In 
the decomposition of X in Eq. (8), all states H for which the energy gap |Mjy — (muo+mu*o)| 
is greater than A^/ {2^£,d*) can be excluded. The effects of such highly virtual states can be 
taken into account indirectly through the scattering length and through other effects on the 
DD* states. If a is large enough, we can choose the ultraviolet energy cutoff A^/ (2/i £)£>*) to 
be smaller than the smallest energy gap. In this case, only the D^D*^ and D*°D^ states in 
the decomposition in Eq. (8) remain and we can set Z = 1. 

The large DD* scattering length could arise from a fortuitous fine-tuning of the pion- 
exchange potential between D^ and D*^ and between D*^ and D^ in the C = + channel. If 
a > 0, there is necessarily a bound state near the D^D*^ threshold, and it is identified with 
the X. Alternatively, the large scattering length could arise from a Feshbach resonance [37], 
which requires the fine-tuning of the potential in a closed channel so that a bound state 
in that channel is close to the D^D*^ threshold. If the closed channel is weakly coupled 
to D^D*^ and D*^D^ in the absence of the fine tuning, then the fine tuning generates a 
large scattering length when the bound state is extremely close to the D^D*^ threshold. 
The natural choice for the closed channel is cc, in which case the bound state in the closed 
channel could be the P-wave charmonium state Xci(2P). Since the decay of the X(3872) 
seems to proceed through J/4> p* and J/ipuj*, one might be tempted to identify the closed 
channel responsible for the Feshbach resonance with J/ipp or J/tpu. However the widths 
of the particles in the closed channel provides a lower bound on the width of the Feshbach 
resonance. The identification of the closed channel with J /ip p or J /ipuj is therefore excluded 
by the fact that the widths of p and lu are larger than the experimental upper bound on the 
width of the AT (3872). 

III. THE DECAY B D^D*°K 

We proceed to apply the separation of the long-distance scale a from the shorter distance 
scales of QCD to the decay process B — ^ D^D*^K. We denote the 4-momenta of the B, D^, 
D*^ , and K by P, g, g*, and A;, respectively. We take the relative 3-momentum q = —q^ in the 
D^D*^ rest frame to be smaller than the separation scale A. The amphtude for B — > D^D*^K 
can be decomposed into three terms corresponding to the three diagrams in Fig. 1. The 
first diagram represents the direct production of D^D*^ by the decay B — > DD*K at short 
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FIG. 1: Feynman diagrams for the decay B D^D*^K, with the short-distance decay amphtudes 
represented by dots and the fong-distance scattering amphtudes represented by blobs. 

distances. The second diagram represents the decay B — >■ D^D*^K at short distances 
followed by the elastic scattering of D^D*^ at long distances. The third diagram represents 
the decay B D*^D^K at short distances followed by the scattering of D*^D^ into D^D*^ 
at long distances. The expression for the amplitude is 

A{B D'^D*'^K] = Abort [5 ^ D''D*''K] 

-I l^,As^ort[B^D'D*'K] 

xD{q + £, moo) D{q, - £, mD*o) A[D^D*' 

The propagators of the virtual D and D* are 

D{p,m) = {p^ - + ie)~^. (14) 

In the loop integrals, there is an implicit ultraviolet cutoff \£\ < A on the 3-momenta of the 
virtual D and D* in the DD* rest frame. 

The long-distance scattering amplitudes in Eqs. (13) are given by the universal expressions 
in Eqs. (9). In the short-distance decay amplitudes in Eq. (13), we can neglect the relative 
3- momentum £ of the D and D*, since it is small compared to all the other momenta in 
the process. The 4-momenta of the D and D* are well-approximated by {rnDo/mx)Q^ and 
(m£)*o/mx)<5'^, where Q = P — k. Lorentz invariance then constrains the short-distance 
decay amplitudes to have the very simple forms 

Ahort[5 ^ D''D*'K] = ci(A) P-e\ (15a) 
Ahort[5 ^ = C2(A) P ■ e\ (15b) 

where e is the polarization 4-vector of the D* . The coefficients ci and C2, which depend on 
the separation scale A, have dimensions of inverse energy. 

In the loop integrals in Eq. (13), the integral over the variable £o can be evaluated by 
applying the residue theorem to the appropriate pole in one of the D meson propagators: 



n*Oi 



(13) 
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FIG. 2: Invariant mass distribution for D^D*^ near thresliold (solid lines) for two values of the 
large scattering length that differ by a factor of 2. The peak in the invariant mass M = mj^o + 
mD*o + /{2ijldd*) occurs at g = 1/a. The crossover from the universal curves to the phase space 
distribution (dashed line) has been modeled by a sudden transition at q = A^. 

where ^(£) is the universal wavefunction of the X given in Eq. (6). The remaining integral 
must be evaluated using the ultraviolet cutoff |£| < A: 

Putting all the ingredients together and keeping only the leading terms for oA » 1 in each 
of the three contributions, the decay amplitude in Eq. (13) reduces to 

A[B D''D"-K] ^ Lw - ^■'A[c.(A) + c.(A)]\ ^ ^. 

\ 'K{l + iaq) J 

Our original expression for the decay amplitude in Eq. (13) simply corresponds to a 
separation of scales, so it is necessarily independent of the arbitrary scale A. However we 
have approximated the long-distance scattering amplitudes by ignoring terms suppressed by 
l/(aA), in which case they reduce to the universal expressions in Eqs. (9). We have also 
approximated the short- distance amplitudes in Eqs. (15) by ignoring terms suppressed by 
i/A. Thus we can expect our expression for the amplitude in Eq. (18) to be independent 



A 



arctan(aA) 



(17) 
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of A only up to terms that are suppressed by powers of 1/A. This is possible only if the 
coefficients Ci(A) scale like 1/A, which implies that the combinations Aci(A) must approach 
ultraviolet fixed points q as A increases. The approach to the fixed points is of course 
ultimately interrupted by the physical scale m^. Neglecting terms that are suppressed by 
powers of 1/A, the decay amplitude reduces to 

A\B ^ D''D*^K] = _^(^l±f^ p . (19) 

The resulting expression for the differential decay rate with respect to the DD* invariant 
mass M is 

j^[B^DD K]-\c, + c,\ 64^5^1^^ 1^^' (20) 

where q is the momentum of the D or D* in the DD* rest frame, 

_ A^/^(M,mgo,mg*o) 
^" 2M ' ^^^^ 

and A(a;, y, z) is the triangle function: 

\{x, y,z)^x^ + y^ + z^ - 2{xY + y'-^" + z'^x^). (22) 

Although we assumed a > in the derivation, our final result in Eq. (20) is valid for either 
sign of a. If a < 0, the only difference in the derivation is that the integral in Eq. (16) 
cannot be interpreted in terms of a universal bound-state wavefunction, since there is no 
bound state. Near the threshold, the DD* invariant mass can be approximated by 

M c:^ rrioo + mo*o + /{2iJ,DD*)- (23) 

The invariant mass distribution in Eq. (20) has a peak at q = l/\a\ with a height that scales 
like \a\ as \a\ increases. It decreases to half the maximum at g = (2 ± V3)/|a|. The full 
width in M at half maximum is 4-\/3/(a^dd*o^). 

The decay rate for B — > D*^D^K also proceeds by the Feynman diagrams in Fig. 1, 
except that the charm mesons in the final state arc D*^ and Z)°. The differential decay 
rate for B — > D*^D^K in the scaling region q <^ is given by exactly the same universal 
expression in Eq. (20). If the large scattering length occurred in the channel with charge 
conjugation C = —, the only difference would be that the factor |ci + in Eq. (20) would 
be replaced by |ci — C2p. We have not specified the charge of the B meson. The fixed-point 
coefficients ci and C2 in Eq. (20) have different values for the decays B~^ —>■ D^D*^K^ and 

The expression for the differential decay rate in Eq. (20) apphes only in the scahng 
region q <C m^r. At larger values of q that are still small compared to the scale — 2mc, 
the resonant terms disappear and the decay amplitude reduces to the short-distance term 
ci(A)P • e* in Eqs. (15a). The corresponding invariant mass distribution dT/dM for q just 
above the scaling region follows the phase space distribution, which is proportional to q in 
the limit g — > 0. The crossover from the resonant distribution proportional to a^q/ {1 + a'^q'^) 
to the phase space distribution proportional to q occurs at a momentum scale that we will 
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FIG. 3: Feynman diagrams for the decay B XK, with the short-distance decay amphtudes 
represented by dots and the long-distance coalesence amplitudes represented by blobs. 



denote by A^. We expect A^^ to be comparable to m^. Just above the crossover region, the 
differential decay rate can be approximated by 

A crude model of the crossover from the phase space distribution in Eq. (24) to the resonant 
distribution in Eq. (20) is a sudden but continuous transition at g = A^, as illustrated in 
Figure 2. This requires 

|c,(A.)| « (25) 

The integral of dV/dM over the region < g < A^^ increases with a. In the limit a —>■ oo, 
it is 3 times larger than the integral of a phase space distribution normalized to the same 
value at g = A^^. 

The Babar collaboration has measured the branching fractions for B D^D*^K and 
B — > D*^D^K using a data sample of about 8 x 10'' BB events [38]. The strongest signal was 
observed in the channel B^ D*^D^K^: 221 ±27 events above the background, but with a 
contamination of about 37 events due to crossfeed from other decay channels. If the invariant 
mass distributions could be measured with resolution much better than m^/(2/i£)£)*) ~ 10 
MeV and if the histograms included enough events, one could actually resolve the resonant 
enhancement near threshold that is illustrated in Figure 2 and determine both a and | Ci + C2 P 
directly from the data. The resolution that would be required may not be out of the question, 
since Babar has presented a histogram of dV/dM for the decay B^ —>■ D*~D*^K^ with 20 
MeV bins [38]. However the region q < m-,^ in which the enhancement is expected to occur 
accounts for only about 0.2% of the available phase space for the decay B —>■ D^D*^K. Even 
with an enhancement in this region by a factor of 3 from a very large scattering length, it may 
be difficult to accumulate enough events in this region to resolve the structure in Figure 2. 



IV. THE DECAY B XK 



We proceed to apply the separation of the long-distance scale a from the shorter distance 
scales of QCD to the decay process B — > XK. The amplitude for B XK can be 
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FIG. 4: Feynman diagrams for the decay B XK: (a) the diagram corresponding to the pro- 

1/2 

duction of the hadron H at short distances, with the blob representing the probabihty factor Z^^ , 
(b)-(c) equivalent diagrams obtained by iterating the bound-state equation to get DD* components 
of the wavefunction. 

decomposed into two terms corresponding to the two diagrams in Fig. 3. The first diagram 
represents the decay B — »• D^D*^K at short distances followed by the coalescence of D^D*^ 
into X at long distances. The second diagram represents the decay B D*^D^K at short 
distances followed by the coalescence of D*^D^ into X at long distances. We denote the 4- 
momenta of the S, X, and K by P, Q, and fc, respectively. The expression for the amplitude 
is 

A[B^XK] = -I I^^A.^ort[B^D'D*'K] 

xD{q + i, rriDo) D{q, - i, m^^.o) A[D^D*^ X] 

xD{q + i, moo) D{q, - i, mD*o) A[D*^D^ X], (26) 

where q'^ = {rajjo /■mx)Q^ and = {m£,*o/mx)Q^ are 4-momenta that add up to the 
4-momentum of the X. 

One might ask why we do not include the diagram in Fig. 4(a), which represents the 
direct production of X through the decay B —>■ HK at short distances, where H is one of 
the hadronic states that appears in the schematic decomposition of the wavefunction of X 
in Eq. (8). Such a short- distance term could be expressed in the form ^short[-B — > HK] 
and can be interpreted as a contribution from a "core component" of the X [36] . The reason 
such a diagram need not be considered is that it is already taken into account through 
the diagrams in Fig. 3. The various components of the wavefunction for the bound state 
X satisfy coupled integral equations. By iterating the integral equations, one can always 
eliminate the component H in terms of a DD* component. Thus the diagram in Fig. 4(a) 
can be expressed in terms of diagrams with explicit DD* states as in Figs. 4(b) and (c). In 
the separation of the amplitude into short- distance parts and long-distance parts represented 
by Fig. 3, the propagators for H in Figs. 4(b) and (c) are absorbed into the short- distance 
decay amplitude. 

The long-distance coalescence amplitudes in Eq. (26) are given by the universal expres- 
sions in Eqs. (10). The short-distance amplitudes for the decays into DD*K are given in 
Eqs. (15). The loop integrals in Eq. (26) can be evaluated as in Eqs. (16) and (17). If we 
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keep only the leading terms for aA S> 1, the amplitude has a factor A[ci(A) + C2(A)]. In 
Section III, we deduced the scaling behavior of the coefficients C\ and C2 as A increases: 
ci Ci/A and ci — > C2/A, where c\ and ci are fixed-point coefficients. The amplitude in 
Eq. (26) therefore reduces to 

A\B^XK\ = (ci + C2)(7r3/XDi)*a)-^/2p.^*_ (27) 
The resulting expression for the decay rate is 

r[g^xjq=ic.+c.r f;ff'r (28) 

This formula applies equally well to the decays XK^ and XK^, with the 

only difference being the values of the fixed-point coefficients Ci and C2. The only sensitivity 
to long distances is through the factor 1/a. We can use the expression for the decay rate 
in Eq. (28) to eliminate the fixed-point coefficients from the expression for the differential 
decay rate in Eq. (20): 

— \B ^ = V[B ^ XK] f^^*''^, . (29) 

dM^ ^ ^ V(l + aV) 



V. ANALYSIS OF S ^ D(*)D(*)i<r BRANCHING FRACTIONS 

A prediction of the branching fraction for B — > XK requires the determination of the 
prefactor |ci + C2P in the expression for the decay rate in Eq. (28). That same prefactor 
appears in the differential decay rate dT/dM in Eq. (20) for B — > D^D*^K in the resonant 
region. Thus measurements of the DD* invariant mass distribution in the resonant region 
could in principle be used to predict the decay rate for B XK. However the resonance 
region q < accounts for only about 0.2% of the available phase space for the decay 
B — s> DD*K. It might therefore be difficult to accumulate enough events to determine 
|ci + C2p directly from the data. There is a crossover from the resonant distribution in 
Eq. (20) to the phase space distribution in Eq. (24) at an unknown momentum scale A^. 
The fraction of the phase space in which dT/dM is described by Eq. (24) should be much 
larger than the 0.2% that corresponds to the resonant region. If one could determine the 
prefactor |ci(A^)p in Eq. (24) from measurements of the DD* invariant mass distribution, 
one could then estimate the desired factor |ci -|- £2^ from the relation in Eq. (25), which is 
based on a crude model for the crossover. The estimate will involve the unknown scale A^, 
which is expected to be comparable to m^r. 

Measurements of dT/dM for the decays B DD*K are not available. The Babar 
collaboration has however measured the branching fractions for decays of B~^ (and B~) and 
of 5° (and S°) into D^*W^*^K, where stands for L>°, D+, L>*°, or D*+ [38]. The 
branching fractions are given in Tables I and II. A substantial fraction of 6 ^ ccs decays 
results in D'^*^D'^*^K final states as predicted in Ref. [39]. The sum of the branching fractions 
is (3.5 ± 0.3 ± 0.5)% for B+ and (4.3 ± 0.3 ± 0.6)% for 5°. An isospin analysis of these 
decays has been carried out [40] . We will use this data to make a rough determination of 
the prefactor |ci(A^)p in Eq. (24). 
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TABLE I: Branching fractions (in %) for — > D^*^D^*^K: measurements from Ref. [38], our 
3-parameter fit, and our 7-parameter fit. 



decay mode 


Br [%] 


3-parameter fit 


7-parameter fit 


D^D+K^ 


0.18 ±0.07 ±0.04 


0.17 


0.18 


D*°D+K^ 


0.41 ±0.15 ±0.08 


0.31 


0.31 




0.52 ±0.10 ±0.07 


0.44 


0.45 




0.78 ±0.23 ±0.14 


0.86 


0.88 




0.19 ±0.03 ±0.03 


0.17 


0.18 


D*^D^K+ 


0.18 ±0.07 ±0.04 


0.31 


0.31 




0.47 ±0.07 ±0.07 


0.44 


0.50 




0.53 ±0.11 ±0.12 


0.86 


0.72 


D-D+K+ 


0.00 ±0.03 ±0.01 





0.00 


D-D*+K+ 


0.02 ±0.02 ±0.01 





0.03 


D*-D+K+ 


0.15 ±0.03 ±0.02 





0.03 


D*-D*+K+ 


0.09 ± 0.04 ± 0.02 





0.13 



TABLE IL Branching fractions (in %) for 5° D^*^D^*^K: measurements from Ref. [38], our 
3-parameter fit, and our 7-parameter fit. 



S° decay mode 


Br [%] 


3-parameter fit 


7-parameter fit 


D-D^K+ 


0.17 ±0.03 ±0.03 


0.16 


0.16 


D-D*^K+ 


0.46 ± 0.07 ± 0.07 


0.41 


0.42 


D*-D^K+ 


0.31 ±0.04 ±0.04 


0.29 


0.29 


D*-D*^K+ 


1.18 ±0.10 ±0.17 


0.79 


0.81 


D-D+K^ 


0.08 ± 0.06 ± 0.03 


0.16 


0.16 


D*-D+K^, D-D*+K^ 


0.65 ±0.12 ±0.10 


0.29 + 0.41 


0.29 + 0.46 


D*-D*+K° 


0.88 ±0.15 ±0.13 


0.79 


0.67 




0.08 ± 0.04 ± 0.02 





0.00 




0.17 ±0.14 ±0.07 


+ 


0.02 + 0.02 




0.33 ±0.21 ±0.14 





0.12 



The most important terms in the effective weak Hamiltonian for b — > ccs decays at a 
renormahzation scale of order nif, is 

Hw = ^^c(,K: (Ci Oi + C2 O2) + h.c, (30) 

where Ci and C2 are Wilson coefficients and Oi and O2 are local four-fermion operators: 

Ci = C7^& s^L^c, (31a) 

O2 = S7^&c7i^c. (31b) 
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We have used the notation 7^ = I'^i^ — 75)- Both operators are products of color-singlet 
currents. We make the simplifying assumption that matrix elements of the operators Oi and 
O2 between the initial-state B and the final-state D^*^D^*^K can be factorized into products 
of matrix elements of currents. For example, the matrix elements for decays into D^D*^K 
and D*^D^K are 

{D^D*^K-\nw\B-) = {GF/V^)V:,Vsc[c,{D'\cr^h\B-){D*''K-\-s^^^cm 

+ C2{K-\s^m-) (i^°5*O|c7^^c|0)) , (32a) 

{D*'D^K-\nw\B-) = (G^/y2)v;;Kc(c'i(/^*°|c7^fe|i?-)(5O;r-|s7^^c|0) 

+ C2 {K-\s-i>lh\B-) (D*°D°|c7^^c|0)) , (32b) 

{D'D*'K'\nw\B') = (Gf/v^)K;K,C2(X°|s7^6|S'^)(L>°L>*°|c7^^c|0), (32c) 
{D*'D'K'\nw\B') = (G^/^/2)T/;K,C2(X°|^7^6|S°)(L>*'^L'°|c7^^c|0). (32d) 

The accuracy of the factorization assumption for this process has been discussed in detail 
in Ref. [41]. 

The terms in Eqs. (32) with coefficient C2 are called "color-suppressed" amphtudes, 
because C2 is suppressed by l/N^ relative to Ci. Only the color-suppressed amplitudes con- 
tribute to the decays B'^ D^*^ D^*^ with Z)*-*-* and D^*^ both charged and to the decays 
B^ D^*)D^*^K^ with D^*^ and D^*^ both neutral. In particular, the only contributions 
to the decays of B^ into D^D*^K^ and D*^D^K° are from the color-suppressed amphtudes. 
As is evident in Tables I and II, the branching fractions for these color-suppressed channels 
are observed to be significantly smaller than those for other decay channels. 

Lorentz invariance can be used to reduce each of the current matrix elements to a linear 
combination of independent tensor structures whose coefficients are form factors. Heavy 
quark symmetry provides constraints between the form factors that can be deduced using 
the covariant representation formalism decribed in Ref. [42]. Matrix elements of operators 
with a heavy quark field Q (or Q) and a Qq meson in the initial (or final) state can be 
expressed in terms of a heavy meson field (or Hy) defined by 

Hv - ^[V^l. + iPvls], (33a) 
= K'^l^ + iPh^] (33b) 

where and Py are operators that annihilate vector and pseudoscalar Qq mesons with 
4- velocity v. We also require the matrix elements of operators with a heavy quark field Q 
and a Qq meson in the final state. They can be expressed in terms of a heavy meson field 
that creates a Qq meson with 4- velocity v: 

Hi = ^ [V:^, - ^n75] , (34) 

The relative phase between the and terms has been deduced by demanding that 
vacuum-to- £)^*^r>(*) matrix elements of operators of the form QTQ have the correct charge 
conjugation properties. 
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We now list the expressions for the matrix elements of the currents that follow from 
heavy-quark symmetry. We denote the velocity 4- vectors of the B, D^*\ and Z)*^*^ by V ^ 
V and respectively. We denote the polarization 4-vectors of the D* and D* by e and e, 
respectively. They satisfy v • e = and v • e — 0. The B-to-D^*'> matrix elements are 

{D{v)\-crLh\B{V)) = ^{w){v + Vr, (35a) 
{D*{v,e)\c'y1b\B{V)) = i^{w) [{1 + v ■ V)e^' - {V ■ e)v^ - ie^{v,V,e)] , (35b) 

where the form factor ^ is a function of w = v -V. We have used the notation e^{p, q, r) = 
e'^^°'^Pj,qari3 and the sign convention = +1. The vacuum-to- L>(*)X matrix elements are 

{D{v)K{k)\s^1c\^) = r)i{K)v^' + r]2{ti)k'', (36a) 
{D*{v,e)K{k)\s-f1c\(ll) = -ir)i{K)e^' -i7]2{K) [{v ■ k)e'' - {k ■ e)v^ + ie^{v,k,e)] , 

(36b) 

where the form factors rji and 772 are functions of k, — v - k. The vacuum-to-D^*)r>^*) matrix 
elements are 

{D{v)D{v)\c^^^c\$) = CK)(^ - vT: (37a) 
{D{v)D*{v, e) |ey^c|0) = iC{w') [{1 - v ■ v)e^ + {v ■ e)v^ + ie^{v, v, e)] , (37b) 
{D*{v, e)D{v)\&f'^c\$) = iC(w') [(1 - v ■ v)e'' + (v ■ e)v^ + ie^'iv, v, e)] , (37c) 
{D*{v, e)D*{v, e)|ey^c|0) = CK) [e • - v)" + {v ■ e)e^ - {v ■ e)e^' - ie^iv - v, e, e)] , 

(37d) 

where the form factor C is a function of w' — v ■ v. The B-to-K matrix elements are 

{K{k)\sjib\B{V)) = ui{k') V^' + u;2{k') k^, (38) 

where the form factors cui and UJ2 are functions of k' — V ■ k. In the current matrix elements 
in Eqs. (35), (36), (37), and (38), the heavy meson states have the standard nonrelativistic 
normalizations. To obtain the standard relativistic normalizations, matrix elements involv- 
ing B, D or D, and D* or D* must be multiplied by m]/^, "O^^^, and m^f , respectively. 

The amplitudes for the decays B — > D^*^D^*'>K at leading-order in A^QCD/f^b and 
^qcd/'itT'c are obtained by inserting the current matrix elements in Eqs. (35), (36), (37), 
and (38) into the factorized expressions for the decay amplitudes, such as those in Eqs. (32). 
For example, the amplitudes for the decays into D^D*^K and D*^D^K are 

A[B- ^ D^D*^K-\ = -iGi{v + V)-e 

-i{G2/mB) [v^ ■ k {v + V) ■ e - ■ {v + V) k ■ e + ie{v + V, v*, k, e)] 

+iG3 [{1 — V ■ V*) V ■ e + {v^ ■ V) V ■ e + ie{v, u*, V, e)] 

+i{G4/mB) [{I - V ■ V:,) k ■ e + {v^ ■ k) V ■ e + ie{v, v*, k, e)] , (39a) 

A[B- ^ D*^D''K-] = iGi[{l + V, ■ V)v e - {vv,)V ■ e - ie{v, v,, V, e)] 

+i{G2/mB) [{l+v^-V)k- e - (v^- k)V ■ e- ie^v^, V, k, e)] 
+iG3 [{1 — V ■ v^) V ■ e + {v^: ■ V) V ■ e — ie{v, v^, V, e)] 
+i{Gi/mB) [{1 — V ■ v.^.) k ■ e + {v^: ■ k) V ■ e — ie{v, v^, k, e)] , (39b) 
^[50 ^ ^o£,*o^O] ^ [(1 - ■ t;.) V ■ e + (t;, ■ \/) t; ■ e + ie{v, v,, V, e)] 

+i{G4/mB) [{I — V • v^.) k ■ e + {v^ ■ k) V ■ e + ie{v, v^, k, e)] , (39c) 
^[^0 ^ 2?*o£,o^O| _ iG3^{i-vv,)V-e+{v,-V)ve- ie{v, v,, V, e)] 

^-i{Gi/mB) [{1—v- v*)k ■ e + {v^- k)v ■ e — ie{v, i;*, k, e)] , (39d) 
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where V, v, and v^, are the velocity 4-vectors of the B, D, and D* and e is the polarization 
4-vector of the D* which satisfies ■ e = 0. We have used the notation e{p,q,r,s) — 
s^^^'^PuqiiTaSp. The four independent dimensionless form factors are 

{mBmDomD*o) ' ^(■y • V) r)i{v^ ■ k), (40a) 

G2{{P-q.f) = {GF/V2)V:,V,,C,{mlmDomD*oy/H{v*-V)rj2{v-k), (40b) 

Gs{{P-kf) = {GF/V2)V:,VscC2{mBmDomD*oy/^C{v-v.)u;,{V-k), (40c) 

G'4((P-A;)') = (G'i./V2)K;KeC2(m|mBomB.o)V2^(^;.^;,)a;2(T/.A;). (40d) 

The amplitudes for the other B — > D^*^D^*^K decays are obtained similarly. Isospin sym- 
metry, in addition to the factorization assumption and heavy quark symmetry, can be used 
to express all 24 decay amplitudes in terms of the four form factors Gi{q^), G2{q^), Gz{q^), 
and Gi{q^). 

Wc proceed to use our expressions for the decay amplitudes to analyze the data from 
the Babar collaboration on the branching fractions for B D^*^D^*^K [38]. For simplicity, 
we approximate the form factors Gi{q^) by constants. We can choose the overall phase so 
that Gi is real- valued. After integrating over the phase space, we obtain expressions for the 
branching fractions that are quadratic in the constants Gi and their complex conjugates. 
The Babar data consists of the 12 branching fractions for B^ given in Table I and the 
10 branching fractions for B^ given in Table II. For each of the data points, we add the 
statistical and systematic errors in quadrature. We then determine the best fits for the 
constants Gi by minimizing the for the 22 data points. 

The decays B+ D^*^D^*^K+ with D^*^ and D^*^ both charged and B° Di*)D(*)K^ 
with D^*^ and D^*"^ both neutral have branching fractions that are significantly smaller than 
other decay channels. The only factorizable contributions to their decay amplitudes come 
from the color-suppressed amplitudes with form factors G^ and G4. Their small branching 
fractions motivates a simplified analysis in which G3 and G4 are set to 0. The only parameters 
that remain are the real constant Gi and the complex constant G2- Thus there are 3 real 
parameters to fit the 22 branching fractions. The parameters that minimize the 

a) 



Gi = 1.9 X IQ-^ (41 
G2 = (-21.2 + 5.5i) X 10"^ (41b) 



The fitted value of G2 is about an order of magnitude larger than that of Gi. The branching 
fractions for this 3-parameter fit are shown in Tables I and II. The per degree of freedom 
is 42.0/19 = 2.2. There are 7 decay modes for which the deviations from the data are 
significantly larger than one standard deviation, including B~^ D*^ . 

We have also carried out a fit that allows nonzero values of the color-suppressed form 
factors G3 and G4. If these form factors are approximated by complex-valued constants, 
there arc 7 real parameters to fit the 22 branching fractions. The parameters that minimize 
the are 

Gi = 1.8 X 10-^ (42a) 
G2 = (-21.6 + 5.0i) X 10-^ (42b) 
G3 = (2.6 + OMi) X 10-^ (42c) 
G4 = (-1.5 - 0.7i) X 10"^ (42d) 
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Note that the values of Gi and G2 are essentially identical to those from the 3-parameter 
fit in Eqs. (41). The branching fractions for this 7-parameter fit are shown in Tables I and 
II. The psr degree of freedom is 29.2/15 = 1.9. There are still 4 decay modes for which 
the deviations from the data are significantly larger than one standard deviation, including 

One could of course improve the fits to the branching fractions by allowing for depen- 
dence of each the form factors Gi, G2, G3, and G4 on the appropriate momentum transfer 
q^. However allowing even for linear dependence on would introduce 8 additional real pa- 
rameters. Such an analysis might be worthwhile if Dalitz plots for the decays were available 
and could also be used in the fits. 

VI. PREDICTIONS FOR B XK DECAYS 

In this section, we use the results of our analysis of the branching fractions for B — > 
D(*W^*^K to estimate the branching fractions for the decays B+ XK+ and S° XK^. 
Our strategy once again is to use that data to provide a rough determination of the prefactor 
k'i(Ajr)P in the differential decay rate dT/dM for B D^D*^K in the region near the DD* 
threshold where the DD* invariant mass distribution follows the phase space distribution 
in Eq. (24). The crossover to the resonant distribution in Eq. (20) occurs at an unknown 
momentum scale A^, which is expected to be comparable to m^r. Given a value for |ci(A^)p, 
we can use the relation in Eq. (25), which follows from a crude model for the crossover, to 
estimate |ci + C2p. This value can then be inserted in Eq. (28) to get an estimate of the 
decay rate for B — > XK. 

We first consider the decay B^ XK^, whose branching fraction should be the same 
as for B~ — > XK~ . The coefficient Ci(A7r) for the decay B~ D^D*^K~ and the corre- 
sponding coefficient C2(A7r) for the decay B~ D*^D^K^ can be deduced by matching the 
amplitudes in Eqs. (39a) and (39b) at the DD* threshold to the expressions in Eqs. (15): 

ci(A^) = C2(A7r) = -iGi/ruB + iG2{mB + + mD'')/m\. (43) 

Using the numerical values for Gi and G2 in either Eqs. (41) or Eqs. (42), the estimate from 
Eq. (25) is 

|ci + C2I ^ 1.6 X 10"^ A^/m^. (44) 

Inserting this into the expression for the decay rate in Eq. (28) and dividing by the measured 
width of the 5+, we obtain 

B.[i,^^Xin.2.TxlO-(A.)^(_|_)"\ (45) 

Our previous analysis in Ref. [34] used the four branching fractions for B^ to decay into 
D^D^K+, D^D*^K+, D*^D°K+, and D*^D*^K+ to fit the constants d and G2. The final 
result was identical to Eq. (45) except that the numerical value of the branching fraction 
for A^ = and Ei, = 0.5 MeV was 2.9 x 10~^. The estimate in Eq. (45) is sensitive 
to the unknown momentum scale Ayr at which the invariant mass distribution crosses over 
from the phase space distribution in Eq. (24) to the resonant distribution in Eq. (20). The 
natural scale for A^ may be rriT^, but we should not be surprised if it differs by a factor of 
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2 or 3. Thus the result in Eq. (45) is only an order-of-magnitude estimate of the branching 
fraction. It can be compared to the product of the branching fractions for — > XK~^ and 
X — > J/ip7i~^'jT~ in Eq. (1). Our estimate is compatible with this measurement if J/ip7i~^Ti~ 
is one of the major decay modes of X. If E^, — 0.5 MeV and if we allow for to differ from 
by a factor of 2, the branching fraction for X — > J/i/jTr^n^ should be greater than lO"-*^. 

We next consider the decay B^ — > XK'^, whose branching fraction should be the same as 
for 5° —>■ XK*^. The amplitudes in Eqs. (39c) and (39d) approach as the DD* approaches 
its threshold. Thus our assumptions of factorization and heavy quark symmetry imply that 
ci(A7r) = C2(A7r) = for this decay. We proceed to consider the size of the coefficients that 
would be expected from the violation of these assumptions. The factorization assumption for 
the B D^*^D^*^K amplitudes can be justified by the large Nc limit. Since we have included 
terms up to 0{1/Nc) in the amplitude, we expect the deviations from the factorization 
assumptions to be 0{1/N^) in the amplitude. Violation of heavy quark symmetry would give 
rise to terms of 0(Aqcd/"^c) in the amplitudes. We expect the largest nonzero contributions 
to the coefficients Ci(A^) and C2{At^) to come from violations of heavy quark symmetry. 

To obtain an estimate of the decay rate for B^ — > XK^, we relax the assumption of 
heavy quark symmetry. Lorentz invariance allows three independent tensor structures in 
the matrix elements {DD*\cy'^c\$) and {D* D\c'y'^c\$) , but heavy quark symmetry requires 
those terms to enter in the particular linear combinations given in Eqs. (37b) and (37c). 
Lorentz invariance implies that only one of the three independent terms can be nonzero at 
the DD* threshold: the term in Eq. (37b) and the term in Eq. (37c). Heavy quark 
symmetry constrains the coefficients of and to be iC{w'){l — v • v), which vanishes 
at the threshold. The constraint of heavy quark symmetry can be relaxed by adding to 
the coefficients of e'* in Eq. (37b) and e'^ in Eq. (37c) the term ixC(l), which is nonzero 
at the threshold. This corresponds to adding the terms ix[G3{V ■ e) + G^/niBik ■ e)] to 
the amplitudes in Eqs. (39c) and (39d). In Table II, the 7-parameter fit gives 0.04 for the 
sum of the two branching fractions for B^ to decay into D^D*°K^ and D^D*^K^, which is 
about one standard deviation below the measured value. The complex parameter x can be 
adjusted so that the sum of the two branching fractions is equal to the central value 0.17 
given in Table II. Using the values of G3 and G4 in Eq. (42), the required values of x form 
a curve that passes through the real values x — —1-9 and x = 5.8 and the imaginary values 
X = ±3.3 If X is allowed to vary over the region in which the sum of the two branching 
fractions is within one standard deviation of the central value, its absolute value has the 
range < |x| < 7.3. 

We proceed to make a quantitative estimate of the decay rate for S° XK^. The 

coefficient Ci(A7r) for the decay B^ — > D^D*^K^ and the corresponding coefficient C2(Ajr) 
for the decay B^ D*^D^K^ can be deduced by matching the amplitudes ix[Gz{y ■ e) + 
G^/niBik ■ e)] to the expressions in Eqs. (15): 

ci(A^) = C2(A^) = ixiGs + G4)/mB. (46) 

If we use the estimate in Eq. (25) to deduce the values of |ci + for both B^ XK° and 
S+ XK'^, the ratio of their branching fractions is 

1^ (47) 

The ratio of the lifetimes of the B^ and B^ is 0.921 ± 0.014. If x is allowed to vary over 
the region < |x| < 7.3, the ratio in Eq. (47) ranges from to 8 x 10~^. We conclude 



Bj:[B^^XK^] 




'IG3 + G4I 


2 


Br[B+ XK+] 


Gi - G2{mB + mo + mD*)/mB\^ 
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that the branching fraction for — > XK'^ is hkely to be suppressed by at least an order of 
magnitude compared to that for XK^ . 

VII. SUMMARY 

If the X(3872) is a loosely-bound S-wave molecule corresponding to a C = + superpo- 
sition of D^D*^ and D*^D^, these charm mesons necessarily have a scattering length that 

is large compared to all other length scales of QCD. The X can be produced through the 
weak decay of the B meson into D^D*^K or D*^D^K at short distances followed by the 
coalescence of the charm mesons at the long-distance scale a. We have analyzed the decay 
B XK and the decays of B into D^D*^K and D*^D^K near the threshold for the charm 
mesons by separating the decay amplitudes into short-distance factors and long-distance 
factors. The long-distance factors are determined by a, while the short-distance factors are 
essentially determined by the amplitudes for B D^D*^K and B D*^D^K at DD* 
invariant masses that are a little above the resonance region, which extends to about 10 
or 20 MeV above the threshold. We obtained a crude determination of the short-distance 
amplitudes by analyzing data from the Babar collaboration on the branching fractions for 
B D^*^D^*^K using a factorization assumption, heavy quark symmetry, and isospin sym- 
metry. 

Our estimate for the branching fraction for S+ — > XK'^ is given in Eq. (45). It scales with 
the binding energy Ei, of X as eI^'^. It also scales as A^, where A^ is an unknown crossover 
momentum scale that is expected to be comparable to m^. If we take E^ = 0.5 MeV and if 
we allow Ayr to vary between and 2m^, our estimate of the branching fraction varies 

from about 7 x 10^^ to about 1 x 10"*^. This range is compatible with the measured product 
of the branching fractions for S+ — > XK'^ and X — > J/'^tt+tt" if Br[X — > J/z/^tt+tt"] is 
greater than about 10~^. 

Our result for the ratio of the branching fractions for B^ XK^ and B^ XK^ is 
given in Eq. (47). It is expressed in terms of parameters Gi, G2, G3, and that appear 
in the amplitudes for B — > D^*^D^*^K. The result is independent of the binding energy E}, 
of the X and also independent of the crossover scale A^. Based on the determination of 
the parameters Gi from our analysis of the Babar data, we concluded that the ratio of the 
branching fractions should be less than about 8 x 10~^. The suppression of 5° — > XK^ can 
be explained partly by the decays of B^ into D^D*^K^ and D*^D^K^ being dominated by 
color-suppressed amplitudes and partly by heavy quark symmetry forcing these amplitudes 
to vanish at the DD* threshold. 

In our analysis of the Babar data on the decays B D^*^ D^*^ K , we made the crude 
assumption that the form factors Gi(g^), G2{q^), Gz{(f )-i and G^^q^) are constants. The 
primary reason for this assumption was that the available experimental information was 
limited to branching fractions for the decays B — > D^*^ D^*^ K . Measurements of Dalitz plot 
distributions and invariant mass distributions for those decays would allow a more rigorous 
analysis that takes into account the g^-dependence of the form factors. This could be used 
to make a more precise prediction of the ratio of the branching fractions for B^ — > XK^ and 
S+ — > XK~^. Measurments of the invariant mass distributions for the decays B — > D^D*^K 
and B D*^D^K would be particularly valuable. They might reveal the enhancement near 
the D^D*^ threshold that would confirm the interpretation of the X as a DD* molecule. 
Even without sufficient data to resolve the peak near the D^D*^ threshold, those invariant 
mass distributions could be used to constrain the parameter A^r in our crude model of the 
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crossover from the resonant distribution to the phase space distribution. This could be used 
to sharpen our estimate of the branching fraction for — > XK^, since the expression in 
Eq. (45) depends quadratically on A^. Measurements of the invariant mass distributions 
would also provide motivation for developing a more accurate model of the crossover. 

The suppression of the decay 5° ^ XK^ compared to B+ XK+ is a nontrivial 
prediction of the interpretation of X(3872) as a DD* molecule. This prediction stands in 
sharp contrast to the observed pattern of exclusive decays of and B~^ into a charmonium 
H plus K. The ratios of the branching fractions for S° — > HK'^ and B'^ — > HK'^ for the 
charmonium states rjc, J/ip, 4^{2S), and Xci(l-P) are 1.33 ± 0.60, 0.85 ± 0.06, 0.91 ± 0.12, 
and 0.59 ± 0.20, respectively. Because charmonium is an isospin singlet and the weak decay 
operators in Eqs. (31) are also isospin singlets, isospin symmetry implies that the ratio of 
the branching fractions for — > HK^ and B~^ — > HK^ should be equal to the ratio of the 
lifetimes r[-B°] and t[-B"'"], which is 0.921 ±0.014. The observed deviations from this lifetime 
ratio are all less than 2 standard deviations. If X were an isosinglet, isospin symmetry would 
imply that the ratio of the branching fractions for B^ XK^ and B^ XK^ should also 
be equal to t[B^]/t[B^]. Thus the observation of suppression of B^ — > XK^ relative to this 
prediction would disfavor any charmonium interpretation and support the interpretation of 
X as a DD* molecule. 

S. Nussinov, who was a coauthor of our previous paper on the decay B^ XK^ [34], 
helped formulate the ideas on which this paper is based. We thank J. Bendich for pointing 
out that a small ratio of the branching fractions for B^ — > XK^ and B^ — > XK^ indicates 
a severe violation of isospin symmetry if X is a charmonium state. This research was 
supported in part by the Department of Energy under grant DE-FG02-91-ER4069. 
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